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Vibrational Frequencies of Semirigid Molecules: 
A General Method and Values for Ethylbenzene ] 

By William J. Taylor 2 and Kenneth S. Pitzer 3 

It is shown that in normal coordinate calculations of the vibrations of molecules it is 
sometimes advantageous to calculate the kinetic-energy matrix, rather than the reciprocal 
kinetic-energy matrix. Explicit formulas are given for the elements of the kinetic-energy 
matrix. Illustrative calculations are given for propane, toluene, and ethylbenzene. A 
semiempirical assignment of the vibration frequencies of ethylbenezene is made on the 
basis of these calculations and the spectroscopic data. 



I. Introduction 

For the purpose of calculating the thermo- 
dynamic functions (particularly the entropy) of 
molecules by statistical methods, it is important 
that the low frequencies of vibration be assigned 
as accurately as possible. When these frequen- 
cies have not been observed in the spectra, the 
assignment must be made either by analogy with 
other molecules, or on the basis of a normal 
coordinate calculation, using force constants 
determined from other molecules. Even when 
the frequencies have been observed, a preliminary 
assignment by analogy or calculation is often 
necessary for the identification of the observed 
frequencies. However, a complete normal coor- 
dinate analysis for molecules with many atoms is 
so time-consuming as to be impracticable in most 
investigations. Therefore, an approximate meth- 
od of calculating the low skeletal frequencies of 
large molecules should be of considerable value in 
connection with thermodynamic calculations. 
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In this paper there is described an approximate 
method that is based on the replacement of the 
actual molecule by a semirigid model (of the same 
dimensions and masses) in which certain groups 
of atoms move as rigid units. In this way the 
number of internal coordinates of the molecule? 
may be reduced sufficiently to make the calcula- 
tion of the low skeletal frequencies practicable. 
The semirigid model should be satisfactory, 
provided that all the frequencies of internal 
motion of the groups assumed to be rigid are 
either considerably higher than the frequencies 
to be calculated, or else do not interact seriously 
with them because of molecular symmetry or 
other factors. In order to reduce the error due 
to interactions, it is desirable that the necessary 
force constants be determined by analogous 
(semirigid model) calculations on other molecules 
in which the frequencies are known. 

The analysis given in this paper shows that 
the kinetic-energy matrix of a semirigid molecule 
is usually more easily calculated than the reciprocal 
kinetic-energy matrix. The reverse is true for a 
complete normal coordinate calculation for the 
molecule. Formulas are given that reduce the 
calculation of the kinetic-energy matrix to a 
consideration of the masses, moments of mass, 
and moments of inertia of the rigid groups. 



An illustrative calculation is given for ethyl- 
benzene by using force constants determined 
from propane and toluene. A semiempirical 
assignment of the vibrational frequencies of 
ethylbenzene is made on the basis of this calcula- 
tion and the spectroscopic data. These calcula- 
tions for ethylbenzene were completed in September 
1945, were and used in the tables of the thermody- 
namic properties of ethylbenzene issued November 
30, 1945, by the American Petroleum Institute 
Research Project 44. 

II. Mathematical Analysis 

1. Kinetic Energy of Semirigid Molecules 

The first step in the normal coordinate analysis 
of a vibrating system is the computation of the 
kinetic-energy matrix, or of its reciprocal. Certain 
general features of the calculation will first be 
considered. Let the system of n particles be 
characterized by 3n generalized coordinates, q, 
six of which refer to over-all translation and rota- 
tion of the system (external coordinates), and the 
remainder of which specify the configuration of 
the system (internal coordinates). The present 
discussion is concerned only with small distortions 
of the system from its equilibrium configuration, 
so that the internal coordinates may pe assumed 
to measure this distortion and to have the value 
zero at the equilibrium configuration. 

The kinetic energy is a quadratic form in the 
coordinate velocities, the #'s, 



2T=qSq. 



U) 



In this equation g is a column matrix, q is a row 
matrix (the transpose of q), and $ is the square 
symmetric kinetic-energy matrix. In general, 
the elements of $ are functions of the q's (not of 
the q's), but because of the restriction to configura- 
tions near the equilibrium configuration, $ is to 
be regarded as a constant matrix here. 

For reasons that will appear presently, it will be 
convenient to assume that the coordinates are 
of three types, and to partition the column 
matrix q accordingly. 



(2) 



The kinetic-energy matrix, $ may be partitioned 
into the corresponding submatrices, and it will 



also be convenient to designate by G the reciprocal 
of $ and to partition G, 
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The concept of a semirigid molecule, with 
which this paper is principally concerned, will 
now be introduced. A semirigid molecule, as 
defined here, is one in which some of the internal 
coordinates, say the group # 3 , remain zero 
throughout the molecular motion, so that g 3 = 0, 
and ^3=0. Thus, if parts of the molecule move as 
rigid units in their equilibrium configurations, the 
internal coordinates of these p'arts will remain 
zero. 

The kinetic-energy matrix for a semirigid mole- 
cule reduces to 



[$11 $12 1 

$12 $22 J 



(4) 



as #3=0. The equivalent expression for the 
kinetic energy is 



2T=q l S ll qi + Hi s nq_2+q.2S 2 2q2- 



(5) 



The form of the potential-energy expression 
will now be considered. As qz=0j these coordi- 
nates do not influence the variable part of the po- 
tential energy (the assumption that q< 6 = is 
equivalent to setting the force constants for these 
coordinates equal to infinity). In addition, there 
will be a group of coordinates, say q x that do not 
occur in the expression for the potential energy, 
although #i occurs in the kinetic energy. This 
group of coordinates of zero frequency includes 
the external coordinates and possibly (as an ap- 
proximation) some very low frequency internal 
coordinates, in particular torsional motions. Final- 
ly, it will be assumed that the potential energy is a 
quadratic form in the remaining coordinates 
(harmonic vibrations) so that, 



2V=q 2 F 22 q 2) 



(6) 



where F 22 is the symmetric potential-energy 
matrix. 

If the Lagrangian equations of a system 
for which the kinetic and potential energies are 
given by eq 5 and 6, respectively, are to be 
integrable, it is necessary to remove by a suitable 
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transformation the cross term between q x and g 2 
in the kinetic energy. The existence of this cross 
term is a consequence of the fact that the initial 
coordinates, q x have been defined in such a way 
that their conjugate momenta, pi=(<)T/dqi), are not 
independent of q 2 . As a result, the momenta p x 
are not conserved, us they should be, because the 
coordinates q { do not enter into the potential- 
energy expression, and are therefore not acted 
upon by any forces. 

As it is desired to find the proper kinetic-energy 
matrix for the (arbitrary) internal coordinates q 2 , 
it follows that the transformation should leave 
these coordinates unchanged. The most general 
linear transformation of this type is 



EH 







a si 



(7) 



where E 22 is the unit matrix with the appropriate 
number of rows and columns. The matrices T n 
and T l2 are subject only to the condition that the 
transformed kinetic-energy matrix have the form 



rs n ] 



(8) 



The equivalent expression for the kinetic energy is 



2T=q[Snq 1 +q 2 S 22 q 2 . 



(9) 



The transformed kinetic-energy matrix is obtained 
by multiplying the original kinetic-energy matrix 
from the right by the reciprocal of the transforma- 
tion matrix, and from the left by the transpose of 
the reciprocal (congruent transformation) . There- 
fore, 



\Sn 1 = [ Tn" 1 

Lo s; 2 \ [-w 
[s 12 s 22 j Lo 





e 22 . 

E 22 , 



Equation 10 is equivalent to the relations 

On= in Oiiin 

0= fn 1 {Su~S 11 Tn l T 12 ) 

S 22 = S 22 — Ti 2 Tu l Si 2 — S r i 2 T 1 r 1 J' 1 2+ 

^12^11 ^11^11 -M2 



(10) 

(11) 

(12) 
(13) 

(14) 



Equations 12 and 13 are equivalent, and yield 

Tn Ti2=Sn S12 (15) 

as the condition that there be no cross term be- 
tween q[ and q 2 in the kinetic energy. 

Equation 15 does not determine T n or T u 
separately, but only the product T^T^, and 
therefore S' n is not uniquely determined. How- 
ever, S22 is unique, as the combination of eq 14 
and 15 yields 



^22 — ^22 ^12*^11 &u 



(16) 



The matrix S 22 will be called the reduced kinetic- 
energy matrix for the internal coordinates (exclud- 
ing those that do not occur in the expression for 
the potential energy). A result equivalent to 
eq 16 has been obtained by Majantz [1] 4 by an 
argument less fundamental than that presented 
here. 

When applied to the external coordinates, the 
l ransformation of eq 7 is equivalent to the applica- 
tion of the conditions that the linear and angular 
momenta of the entire molecule are independent 
of the q/s. When applied to an internal tor- 
sional motion, the transformation is equivalent to 
the condition, previously stated by Sayvetz [2, 3], 
that the momentum conjugate to the torsion 
should he independent of the (remaining) ^'s. 

An important property of transformations of 
the type of eq 7 may now be stated. The proof 
is given in appendix 1 . The coordinates q 1 may 
be divided into two or more groups in various 
ways, and the transformation applied consecu- 
tively to each group. The over-all transforma- 
tions from q x to q[ will not, in general, be identical. 
However, the theorem states that all possible 
transformations of this type lead to the same 
reduced matrix S' 22 . The importance of this 
theorem is that in practice it is convenient to 
apply the transformation consecutively to trans- 
lation, rotation, and finally, in some cases, to 
certain internal coordinates. 

The Lagrangian equations for the system may 
now be written as 



p'i=Sn([i=0 

S 22 q 2 + F 22 q2=0 



(17) 
(18) 



4 Figures in brackets indicate the literature references at the end of this 
paper. 
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The first equation shows that the momenta p[ are 
conserved. The second equation may also be 
written 



rit ten 
where 



§2+^22^2222 = 0, 

G22=(S22)~ 1 



(19) 

(20) 



is the reduced reciprocal kinetic-energy matrix for 
the internal coordinates q 2 . In eq 16, S 22 is de- 
fined in terms of the submatrices of S. It is 
proved in appendix 2 that 67 22 may be expressed 
in terms of the submatrices of G, eq 3, by the 
equation 



G 2 2 — ($22) — G 22 — ^23 G33 1 G% 



(21) 



This result has been obtained previously by 
Wilson [4], although not explicity for the case in 
which coordinates of the types q x and q z are present 
simultaneously. 

As is well known, eq 18 and 19 lead to the 
following secular equations for the determination 
of the normal frequencies of vibration of the 
coordinates q 2 : 



or, 



\g; 2 f 22 -\e 22 \=o, 



(22) 
(23) 



wmere \=4t 2 v 2 (v = normal frequency). In the 
present method, in which £ 22 rather than 6? 22 is 
computed, it is convenient to convert eq 22 to the 
secular equation of a single matrix, in the form 



\S22F2i 



x -1 ^ 



0. 



(24) 



In practice, F 22 is diagonal, or nearly so, and 
therefore easily reciprocated. The roots of eq 
24 are, of course, the reciprocals of the X's. 

It is possible now to state clearly the relative 
advantages of calculating £ 22 and 6? 22 (see eq 3, 16, 
and 21). In calculating S22, the coordinates q 3 are 
ignored, but all the elements for q x and q 2 must be 
computed. Conversely, in calculating c? 22 , the 
coordinates q x are ignored, but all the elements for 
q 2 and g 3 must be computed. It should also be 
noted that even with the simplifications described 
later in this paper, the elements of S are more 
difficult to compute than those of 6? 5 (this state- 

5 Wilson [4] has described a now well-known method for the calculation of 
G. A detailed exposition of the method has been given by Meister and 
Cleveland [5j. A method similar to Wilson's has been developed indepen- 
dently by the Russian workers M. Eliashevich, B. Stepanov, and M. Wolk- 
enstein. References to the papers of these workers are not given here, as we 
are not directly concerned with the calculation of G. 



merit does not apply to the elements for the 
external coordinates, but these are just the ele- 
ments of G that may be ignored). From these 
considerations it is obvious that the calculation 
of G 2 2 is to be preferred when all or most of the 
internal coordinates are included in the calculation 
(so that the number of coordinates of type q s is 
zero or small). On the other hand, S& is often to 
be preferred in calculations by the semirigid mole- 
cule approximation as the number of coordinates 
of the type g 3 is then relatively large, and the 
number of coordinates remaining in the calculation 
(type q 2 ) is relatively small. 

If the molecule has the symmetry of one of the 
point groups, the internal coordinates fall into sets 
equivalent under the group. The coordinates in 
an equivalent set may be transformed into each 
other by the operations (rotations and reflections) 
of the group. In this case, just as is the case when 
working with the reciprocal kinetic-energy mtarix 
[4, 5], it is desirable to construct new coordinates, 
symmetry coordinates, which are linear combina- 
tions of equivalent internal coordinates, and which 
reduce the group. That is, the symmetry co- 
ordinates are chosen in such a way that each trans- 
forms according to one of the irreducible repre- 
sentations of the group. The advantage of the 
use of symmetry coordinates is, that there are no 
interaction elements in the kinetic- or potential- 
energy matrices between coordinates which trans- 
form according to different irreducible representa- 
tions. The matrices therefore factor into sub- 
matrices (the number of identical submatrices for 
each representation is equal to the dimension, or 
degeneracy, of the representation), and the 
secular equation factors into the corresponding 
equations of lower degree [6, 7]. 

It will now be shown that the symmetry 
coordinates may be constructed either before or 
after applying the transformation of eq 7, provided 
that none of the symmetry coordinates involve 
linear combinations of coordinates of types q x and 
q 2 . With this restriction the symmetry coor- 
dinates may be represented by 



[?•■ lm 



(25) 



The matrices U n and U 22 are usually, although not 
necessarily, orthogonal. The transformed kinetic 
energy matrix is 
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[tin 1 o If s n ^"ir^r 1 o 
Lo u^JLSu s. 22 Jlo c/ 22 -'J 

[Uu'S^Un 1 u 1 r i s 12 u 22 - 1 ~] (2 

The reduced kinetic-energy matrix is then, 
fr -iQ TT - 1 

(_y 22 ^22 ^ 22 

l\l l S 13 Uu » ( % 'Sntftr " ' J7u '&• E7* 1 

= U22 ^22^22 ^22 ^12^11 ^12^22 

= L 22 (^22 ^12^11 ^12.)^22 > 



(27) 



which is identical with the result obtained by 
transforming the original matrix aS 22 of eq 16 by 
£7 22 . It will also be observed that the matrix 
S22 is unaffected by symmetry combinations 
formed with the coordinates q± (such, for example, 
as a rotation of axes). 

When U 2 2 is orthogonal the transformed kinetic- 
energy matrix is U^^Un, but for nonorthogonal 
U22 the general form U 2 2 l S 22 U2 2 ~ l must be used. 
An illustration of the use of a nonorthogonal 
transformation will occur in the calculations at 
the end of this paper. 

2. Calculation of the Kinetic-Energy Matrix 
(a) Equations for Generalized Coordinates 

Explicit formulas for the calculation of the ele- 
ments of the kinetic-energy matrix for generalized 
internal coordinates are presented in this section. 
Denote by m* the mass of the ith atom, and by 
R* the position vector of the ith atom with respect 
to a point stationary in an inertial system. Then 
the kinetic energy of the system of particles is 



2T=]TXR r R„ 



(28) 



where R* is the vector velocity of the it\\ atom. 

It will be assumed here, as in the preceding 
section, that the positions of the n particles are 
specified by 3/i generalized coordinates q, some of 
which, as before, remain zero for a semirigid 
molecule. 6 We may then define a set of vectors, 

**,=(dR</aa«). (29) 



8 In this section specific coordinates will be indicated by letter subscripts, as 
g m and g n . 



Thuss^iis a vector pointing in the direction in 
which the ^th atom moves when q m increases from 
zero (for the equilibrium configuration) to a small 
positive value, all other g's remaining zero, and the 
magnitude of s mi is equal to the displacement of 
the ith atom per unit change in q m . Then, 



R?' — Z-fimiQmy 



(30) 



and 



i \ m / \ 11 


4n\ 




m,n \ i ) 




(31) 


But from eq 1, 






2 T= *s,S m . n q m q_ n , 

m,n 




(32) 



where S mn is the element in the mth row and nth 
column of the symmetric kinetic-energy matrix 
S. It follows from eq 31 and 32 that 



*^ mm x j Wl i S m i • S m j , 

i 

& mn — ^ n m ~~ / i "» /S m i ' S^ j 



(33) 

(34) 



In order to proceed further, define an orthogonal 
system of x, y, z, -axes in the molecule, and denote 
the unit vectors, directed positively along the 
axes, by x, y, and z. Let the system be right- 
handed, so that xXy = z. The position of the 
origin of the system, and the directions of the 
unit vectors, must be uniquely determined with 
respect to the molecule for any arbitrary (small) 
displacements of the internal coordinates, but the 
choice of axes is otherwise arbitrary. In parti- 
cular, it is not necessary that the origin be at the 
center of mass of the molecule. 

The coordinates of the ith atom are x u y u and 
2i. For a given choice of axes these coordinates 
are completely determined when the values of the 
internal coordinates are given. The position 
vector of the iih atom with respect to the origin is 



Ti^xXi+yyi+iZi. 



(35) 



Three of the generalized coordinates will be 
taken to be the translations of the origin of the 
axes in the (instantaneous) directions of the x-, y-, 
and z-vectors. The r-vectors of eq 29 are identi- 
cal for each atom, for each translation, and are 
simply the unit vectors 



s x — x, s^ — y, s 2 



(36) 
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Three more coordinates will be taken to be the 
three rotations, by angles <j> x <t> v and <j> z about the 
x- ) y-, and 2-axes. The S-vectorsfor these are easily 
shown to be the vector products, 

For the remaining coordinates, the (3n— 6) in- 
ternal coordinates, 



s mi =(dr</dg m ) 



(38) 



where, 



(dr i /d(2 w )=x(b;r ? /5(/ m ) + y(d^/,-/dg w ) + z(d0 i /dg w ). 

(39) 

Explicit formulas for the elements of the 
kinetic-energy matrix may now be obtained from 
eq 33 and 34, combined with eq 36 to 39. The 



reductions of the scalar products of the s- vectors 
are made by the standard methods of vector 
analysis. The reduction process of eq 16 may 
then be applied explicitly to the translations. In 
this reduction, S n in eq 4 and 16 represents the 
submatrix for translation, and is a 3X3 diag- 
onal matrix, with each diagonal element 
equal to the total mass of the molecule, M. 
$22 represents the submatrix for rotations and 
internal coordinates, and S i2 the submatrix 
for the interaction of these coordinates with 
translation. The elements of S l2 are simply 
the X-, y-, and E-components of linear momentum 
of the rotations (about the x, y, 2-axes) and 
internal coordinates. 

The reduced matrix for the rotations and in- 
ternal coordinates, corresponding to S 22 in eq 16, 
which is obtained in this way may be written, 



<f>y 



q n 



4>x 


T x 


J xy 


-hx ■ 


. H xm . . 


. H 


<t>v 


-• xy 


h 


J yz 


£1 ym . . 


. H 


<t>z 


~hx 


-/» 


h • 


H zm . . 


. H. 



H x 



H x 



Hy 



Hy, 



H„ 



H„ 



K m 



K nn 



K m 



K n 



(40) 



The elements for rotation are the moments and 
products of inertia of the molecule about axes 
parallel to the j,7/,2-axes, but passing through 
the center of mass of the molecule, 

I x =^mM + z$) -M(yZ+ zl) 
Iy^m^zl+xD-Mizl+xl) 
It^mM+yV-Mixh+yh) 

I xy = ^m iX t y x — Mx M y M 

Ivz=5fmiyiZi—My M z M 

I zx =^m l z i x i —Mz M XM 

The total mass of the molecule has been denoted 
by M, and the coordinates of the center of mass 
of the molecule, for the equilibrium configuration, 



by x M , y M , and z M . The first term in each equa- 
tion of 41 is the moment or product of inertia 
about the arbitrary origin; the second term arises 
from the reduction process of eq 16, and vanishes 
if the origin coincides with the (equilibrium) 
position of the center of mass. 

The elements for the interaction of rotation and 
internal coordinates may be written 



(41) 



= z-H» 



(42) 



tt-xm X , ri m , tlym — y*xl m , £1 z\ 

where H m is the vector 

Hm^^miTiX (i>r i /iki m )—r M X^;m. i (dr i /d(i m ) (43) 

i i 

The vector from the origin of the axes to the center 
of mass of the molecule has been represented by 

r M =M- l J2m i r i =xx M +yy M +'zzM. (44) 
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If the origin of the x,y,z-&xes coincides with the 
(equilibrium) position of the center of mass of the 
molecule, the vector r M , and the second term in 
eq 43, vanish. The second term arises from the 
reduction process of eq 16. The first term in eq 
43 represents the vector angular momentum of the 
motion g m about the origin, while H w represents 
the vector angular momentum of q m about the 
(equilibrium) center of mass. The quantities 
II xm , Hym, and H zm are therefore the components of 
angular momentum of q m about axes through the 
center of mass parallel to the x, y, £-axes. 
For internal coordinates, 

^m=^m,(dr,/5 2 m) 2 -M~ i r2:m i (dr,/d2 m )l 2 (45) 

and for the interaction of two internal coordinates, 

K mn =K nm J2m i (dr i /dq m )'(dr i /dq n ) 

i 

-M-T^]m,(5r,/d(Z m )l.r^m,(dr,/dg re )J (46) 

The second term in eq 45 may also be written 
M(dr M /dq m ) 2 , and therefore vanishes if the position 
of the center of mass of the molecule is independent 
of q m . Similarly, the second term in eq 46 is 
M(dr M /dq m )-(dr M /dq n ), and vanishes if the center of 
mass is independent of either q m or q n (or both). 
However, the vanishing of the (equilibrium) 
vector x M is not, of itself, sufficient to cause these 
terms to vanish. The second terms in eq 45 and 
46 represent the kinetic energy of translation of 
the center of mass of the molecule, and arise in 
the reduction process of eq 16. The first terms in 
these equations represent the kinetic energy of 
the internal coordinates computed with respect 
to the arbitrary origin. Subtraction of the second 
term from the first yields K mm and K mn , which 
represent the kinetic energy of the internal co- 
ordinates computed w T ith respect to the instan- 
taneous position of the center of mass. 

The set of atoms associated with the internal 
coordinate q m will be defined to include those 
atoms which move with respect to the x,y,z-eLxes 
when the coordinate q m varies. This set is not to 
be regarded as fundamental, because the atoms 
included in it depend, in general, on the definition 
of the axes. The vector (dri/dq m ) vanishes by 
definition for all atoms not in the set associated 
with q m . As a result, the summations in eq 43 
and 45 need extend only over the set of atoms asso- 
ciated with q m . In eq 46 the two summations in 



brackets extend over the sets of atoms associated 
with q m and q n , respectively, whereas the first 
summation in eq 46 extends only over those atoms 
included in both of these sets. 

The final reduced kinetic-energy matrix is ob- 
tained from the matrix 40 by applying the reduc- 
tion process of eq 16 to the rotations, and to any 
internal coordinates that do not appear in the 
potential-energy expression. The elements of the 
matrix of eq 40 are not independent, in general, 
of the choice of the x,y,z-&xes, but the reduced 
matrix is independent of the axes after the over- 
all rotations have been removed. The reduced 
matrix is most easily obtained numerically, by 
applying the reduction process of eq 10 to one 
row and column of the matrix at a time (so that 
in each step S n in eq 16 is a matrix with a single 
element). This procedure makes it unnecessary 
to reciprocate directly the submatrix for the co- 
ordinates that are to be removed. 

(b) Equations for semirigid molecules 

The equations for the calculation of the kinetic- 
energy matrix for generalized internal coordinates 
have been given in the preceding section. These 
equations may now be applied to the two types 
of internal coordinates which will be found use- 
ful in specifying the configuration of a semirigid 
molecule. These arc coordinates that represent 
either a translation or a rotation of their associated 
sets of atoms with respect to the x,y,z-axe8. In- 
ternal translational and rotational coordinates will 
be represented by q t and q r , respectively, where 
q t is the distance of translation, and q r is the angle 
of rotation, both measured from the equilibrium 
configuration. 

Let ix be a unit vector in the direction of transla- 
tion for the coordinate q t , with components 
n x , Liy, and fi z , and let £ be a unit vector directed 
along the axis of rotation for the coordinate q T) 
with components £ x , £ y , and £ 2 . Let p t be the posi- 
tion vector of the iih atom associated with q r , 
with respect to an arbitrary point on the axis of 
rotation for q r . In terms of these vectors eq 38 
becomes, 

(dr,/d<z,)=M, (47) 

and 

(&r l /&ff r )=jXpi l (48) 

for translational and rotational coordinates, re- 
spectively. Equations 47 and 48 may now be 
substituted in eq 42 to 46, to obtain the elements 
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of the kinetic-energy matrix, eq 40, in vector form. 
These vector expressions may then be reduced by 
the usual methods to more useful scalar expressions 
involving the masses, moments of mass, and 
moments of inertia of the associated sets of atoms. 
The details of these reductions will be omitted and 
only the final formulas stated. 

Some additional definitions will be necessary. 
The total mass of the set of atoms associated with 
a coordinate q m (m=t or r) will be denoted by 
M m , and the coordinates of the center of mass of 
the set of atoms by x Mm , y Mm , and z Mm . The 
corresponding quantities for the entire molecule 
are M, and x M , y M) and z M . 

For each rotational coordinate, q r , define a set 
of orthogonal right-handed £, 77, f-axes, with the 
origin at a point x or , y or , z or , on the axis of rotation, 
but otherwise aribtrary, and the £-axis coinciding 
with the axis of rotation. The direction cosines 
of the £, 77, f-axes with respect to the x, y, 2-axes, 
written in matrix form, are 



% 



u 



Vx 


k: 


Vy 


r, 


Vz 


rJ 



(49) 



The determinant of this matrix will be + 1 if 
both the x, y, z-, and J, 77, f-axes are right-handed. 
The following moments and products of inertia 
of the set of atoms associated with q r will be 
required, 

J«=2>«(i?+tf) 



i 



(50) 



1% is the moment of inertia about the axis of rota- 
tion. The products of inertia, 7^ and 7^, vanish 
if (1) the £-axis is a twofold, or higher, axis of 
symmetry for the associated set of atoms; (2) the 
77, f -plane is a plane of symmetry; or (3) if a three- 
fold, or higher, axis of symmetry lies in the 77, 
f -plane, and intersects the £-axis. 7$, (but not 
necessarily 7^) vanishes if (1) the f,£-plane is a 
plane of symmetry or (2) the 77-axis is a twofold 
axis of symmetry. 

Let d r be the perpendicular distance from the 
center of mass of the set of atoms associated with 



q r to the axis of rotation, or £-axis, and let X^, \ y , 
and X 2 be the direction cosines, with respect to 
the x,y,z-&xes, if the line along which this center 
of mass moves when q r increases from zero to a 
small positive value. 
That is, 



and 



dr=(i7i+fi), 

\ = \VMrty~ ^MrVy) /d r 
X 2 = (VMrtz — ^MrV 2 )l(lr, 



(51) 



(52) 



where, 



VMr=M- r 1 ^2m iVi , Ur=M~ r 1 ^2m, i ti. (53) 



The sums in eq 50 and 53 extend over the atoms 
associated with q r , Atoms on the axis of rotation 
are not associated with q r , according to the 
definition given in the preceding section. It may 
be more convenient to calculate these quantities 
as follows. Define a x , a Vj and a 2 to be equal to 
their respective cofactors in the determinant, 



OL X 
\X]\fr %or) 



vJmt V or) 



& 

{2 Mr Zor) 



(54) 



Then 



and 



d r =(a* + a£-f-a?), 



(55) 



*x= (ot x /d r ) , \y=(a !/ /d T ) } \ z =-- (ajd r ) . (56) 

The equations for the elements of the kinetic- 
energy matrix, eq 40, may now be stated. For a 
translational coordinate, q t , H xt , H yh and H zt are 
equal to their respective cofactors in the deter- 
minant 



(57) 



H xt 


H„ 


H lt 


M tf i x 


M„ji V 


M tl x, 


\%Mt %MJ 


\\)Mt — Vm) 


(Zmi—Zm) 
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[f the translation arises from the stretching of a 
bond, n x . Hy, and ix z are the direction cosines of 
the bond. 



For a rotational coordinate, q r , 



H vr =m 






(58) 



where!!™,!!™, and H^ aredefmed by the equations 
H £ = HJi - vxhv - tJn 
FS?=^/ l - % J^-f,/ rr , (59) 

* andflf^, H ( yjl, and H„ are equal to their respective 
cofactors in the determinant 



H {2) 

±J - xr 

M r d T \ x 



M4r\ 



J J (2) 

M r d T \ z 



(60) 



In many cases, 7 €l? =/^=0, and eq 59 reduces to 

ffS?=&/«, H5={A H2?=fc/|. (61) 

Before stating the formulas for the remaining 
elements, iT mw , and i^ mn , of eq 40, it will be neces- 
sary to differentiate between connected and un- 
connected coordinates. If the sets of atoms asso- 
ciated with two coordinates q m and q n contain no 
atoms in common, the two coordinates will be 
said to be unconnected — otherwise they are con- 
nected. Whether two coordinates are connected 
or unconnected will be dependent, in general, on 
the way in which the x,y,z-&xes are defined. The 
only types of connected coordinates that will be 
considered are those for which the set of atoms 
for one coordinate, say q n , is completely contained 
within the set of atoms for the other coordinate, 
q m . A special case is that in which the sets of 
atoms for q m and q n are identical. If q m and q n 
are unconnected, this fact will be indicated ex- 
plicitly by writing K mn {=K nm ) as K mtn (=K n , m ). 
If q m and q n are connected, and the set of atoms 
for q n is completely contained within the set of 
atoms for q m , then K mn (=K nm ) will be written 
K m{n) (=K (n)m ). If the set of atoms for q m and 
q n are identical, K m{n) =K [m)n . Subscripts t and r 
will again be used to indicate translational and 



rotational coordinates, and a second coordinate 
of a given type will be distinguished by the use of 
primes. The equations for the several possible 
forms of K mm and K mn may now be written as 

K lt = M t {l-M-'M t ) (62) 

K ttl >= -M-'MMt'iw'x+w'y+w'z) (63) 

K nn = M/(l^M-W t )(^ x + ^, + ^^) (64) 

K t , r =-M- l MMAtf (65) 

K tir) =M r (l-M-'M t )A tr (66) 

K rit) = M t A r(t) -M-*M t M r & )r (67) 

K rr =It-M- l M 2 r d 2 r (68) 

K r , T >= -M- l M r M' r d r d' r (\ x \' ' x +\,y ' V +W f ',) (69) 

-~M- l M T M' r d r d' r (\ x \' x + XX+ KK) (70) 

The symbols A tT , A T ( t) , and A r ( r '), represent the 
following determinants: 



A tr = d r (n x \ X -\- Hy\y+ jJL g \ Z ) 









Vz 



[Xmt Xor) \Vmt Vor) (<Mr Z or ) 



(71) 



&Tit) = 



fJ-x 



fXy 



V>z 

t 



\$Mt %or) (]hlt. Vor) \%Mt Zor) 



(72) 



A r{r 



(r')" 



K 



K 



K 



(XoT—X or ) (y or — Vor) (z' or —Z or ) 



(73) 



In eq 70, #<»', Hg>', and HJ»', are the quantities 
calculated from eq 59 (with r—r f ). When eq 61 
applies, the first term in parentheses in eq 70 
reduces to, 

(u:+^;+«:)4 (74) 

If the center of mass of the set of atoms asso- 
ciated with a rotation, q r , lies on the axis of rota- 
tion, so that d r =0 y q r will be called a balanced 
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rotation. For a balanced rotation, H™=H™= 
H { ^ = 0, in eq 58, and there are obvious simplifica- 
tions in eq 65 to 70, as A, r =0 (in addition to 
d r =0). A balanced rotation may be called acci- 
dentally balanced if the £-axis is not an axis of 
symmetry of the associated set of atoms. Such 
cases are very unlikely to occur. As a rule, the 
£-axis for a balanced rotation will be a twofold or 
higher axis of symmetry. In this case, 1^=1^ = 0, 
and eq 59 is replaced by eq 61. However, it is not 
necessary that a rotation be balanced for /$, and 
7f$ to vanish. 

III. Application to Propane, Toluene, and 
Ethylbenzene 

1. Kinetic-Energy Matrices 

The application of this method to certain vibra- 
tions in ethylbenzene will be described in this 
section. In order to obtain the necessary force 
constants for ethylbenzene, it also has been neces- 
sary to set up the kinetic-energy matrices for 
propane and toluene. In the semirigid model of 
ethylbenzene for which the calculation has been 
carried out, it is assumed that the benzene ring 
and the attached hydrogens move as a rigid body, 
and also that the methylene and methyl groups 
move as rigid bodies. That is, the force constants 
for all deformations of the benzene ring and 
attached hydrogens and for deformation of the 
H — C — H angles in the methylene and methyl 
groups have been set equal to infinity. The force 
constants for the stretching of the two C — C bonds 
in the ethyl group have also been assumed infinite. 7 
In order that the force constants obtained from the 
observed frequencies of propane and toluene 
should be on a comparable basis, corresponding 
assumptions have been made for these molecules. 
The coordinates remaining in the problem are 
therefore the rocking and twisting motions of 
these rigid groups. 

The total number of internal coordinators in 
ethylbenzene is 48; of these, 8 are retained in the 
present calculation, and 2 are assumed to have 
zero frequency. In calculating the kinetic-energy 
matrix by the present method it is necessary to con- 
sider only these 10 internal coordinates, plus the 



3 over-all rotations, and then to remove the 5 co- 
ordinates of zero frequency, as discussed in 
section 11,1. In order to obtain the reciprocal 
kinetic-energy matrix for the 8 internal coordi- 
nates, it would be necessary to set the matrix up 
for the 46 internal coordinates of nonzero fre- 
quency, and then to remove 38 coordinates, as 
discussed in section 11,1. 

Before enumerating the coordinates, the sym- 
metry properties of the molecules will be reviewed 
briefly. The propane molecule has the symmetry 
O2V, with the nondegenerate representations A iy 
A 2 , B u and B 2 . Representations A and B are 
symmetric and antisymmetric, respectively, to 
rotation by 180 degrees about the twofold axis. 
Representations with subscripts 1 and 2 are sym- 
metric and antisymmetric, respectively, to reflec- 
tion in the symmetry plane. The toluene and 
ethylbenzene molecules have the symmetry d h , 
with nondegenerate representations A f and A" , 
which are symmetric and antisymmetric, respec- 
tively, to reflection in a plane perpendicular to the 
plane of the benzene ring. 8 

The notation adopted for the internal coordi- 
nates is as follows. In-plane and out-of-plane 
(with respect to the symmetry plane) rocking 
motions are denoted by <r< and ic i} respectively. 
Internal rotational motions are denoted by </>*. 
The subscript i indicates the group, as follows: 
Propane, methyl (1 and 2), and methylene (3); 
ethylbenzene, methyl (1), ethyl (2), and methylene 
(3); and toluene, methyl (1), and methyl group 
(2). 9 On this basis methylene twisting is repre- 
sented by 3 . The C — C — C angle bending (in 
propane and ethylbenzene) is denoted by a. The 
coordinates represent in each case the angular 
deformation from the equilibrium configuration. 
In addition the over-all rotations about the 
x,y,z-&xes will be denoted by <j> x , <t>y, and <j> z , 
respectively. 



7 As no bond stretchings are included in the present calculations, they will 
not illustrate the use of the formulas for translational internal coordinates 
given in section II, 2b. 



8 The equilibrium configuration assumed for the ethylbenzene moleclue 
is that in which the plane of the benzene ring is perpendicular to the plane 
of the C—C — C angle in the ethyl group. The angle of rotation of the methyl 
groups in these molecules does not affect the quantities entering into the 
calculations. 

9 The axes for methyl and methylene rocking pass through the carbon atoms 
of the CH3 and CH2 groups, respectively, so that these are motions of the 
hydrogen atoms only. The axis for ethyl rocking passes through the phenyl 
carbon atom, sc that the entire C2Hj group moves as a unit. The motion in 
toluene of the entire CH3 group as a unit about an axis through the phenyl 
carbon atom has been called "methyl group" rocking to distinguish it from 
methyl rocking. Methyl rocking may conveniently be denoted by CH3 
rocking, and methyl group rocking by — CH3 rocking. 
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Iii the case of propane it is necessary to con- 
struct the following symmetry coordinates: 



<7 + =(l/V2)(<r l +<r 2 ) 

<r=(l/V2)(«i-«ri) 

7T+=0/V2)(7r 1 + 7r 2 ) 

T-=(1/J2K*1-TJ 

* + =(l/V2)(*i+*») 
*-=(l/V2)foi-fc) 



Table 1. — Notation and symmetry for coordinates 



(75) 



These transformations are orthogonal, with the 
matrix, 



(1/V2) 



1 



1 
-1 



(76) 



\ 



i 



1 



From eq 26 there are obtained the relations 

(cr+(7 + ) = (l/2)[((r 1 ,(T 1 ) + (<r 2 ,(r 2 )]+((r 1 ,(r 2 ) 

(c- J a-) = (ll2)[(a u o' 1 )+(a 2 ,a 2 )]-(a u a. 2 ) 
(cr + ,<r-) = (l/2)[((r 1 ,(r 1 )-((r 2 ,(r 2 )] = 
(o' + ,(z) = (l/V2)[(<r 1 ,2) + ((T 2 ,g)] 

(o--,q) = (lN2)[(<T l ,< 1 )-(o- 2 ,q)] 



where (o- + ,<r + ) is the diagonal element in the 
kinetic-energy matrix for o- + , (<j + , <r~) is the element 
connecting a + and o-" (equal to (a~, o- + )), etc. 
# represents any coordinate other than n, <r 2 , 0" + , 
or a - . Similar relations hold for the w and </> co- 
ordinates. The elements for the internal coordi- 
nates are first calculated from the equations given 
in section II, 2b, and the elements for the sym- 
metry coordinates are then calculated from eq 76. 
The notation used for the coordinates, their 
symmetries or representations, and their conven- 
tional descriptions are summarized in table 1. 
Interaction terms in the kinetic- and potential- 
energy matrices between coordinates belonging to 
different representations vanish identically. 



Symmetry 


Coordinate ° 


Description 


PROPANE (C 2 .) 


Ax 


a 
a 


CH3 rocking 
C-C-C bending 


A% 




Over-all rotation 
CH 3 rotation 


7T + 


CH2 twisting 
CH3 rocking 


B x 


<Py 


Over-all rotation 


<T + 


CH3 rocking 
CH2 rocking 


B 2 


<px 


Over-all rotation 
CH3 rotation 


1T~ 
7T3 


CH3 rocking 
CH2 rocking 


ETHYLBENZENE (C u ) 


A' 


Vv 


Over-all rotation 


a 


CH3 rocking 
C2H5 rocking 
CH2 rocking 
C-C-C bending 


A" 


<Px 
<Pz 

<P2 


Over-all rotation 
Over-all rotation 
CH3 rotation 
C2H5 rotation 


7T2 
7T3 


CU2 twisting 
CH3 rocking 
'ii rocking 
CH2 rocking 


TOLUENE (C 1A ) 


A' 


<Py 


Over-all rotation 


<T2 


CH3 rocking 
-CH3 rocking e 


A" 


<Px 


Over-all rotation 
Over-all rotation 
CH3 rotation 


7T1 
7T2 


CH3 rocking 
-CH3 rocking 



a The coordinates above the dashed line in each representa tion were 
removed from the kinetic-energy matrices (see text). 
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The dimensions used in the calculations are as 
follows. The benzene ring was assumed to be a 
plane hexagon, with the length of the C — C bonds 
1.39 angstrom units. The other bond lengths, in 
angstrom units, were 1.54 for the C — C bonds in 
propane and in the alkyl groups (including the 
bond joining the alkyl group to the ring [8]), 1.09 
for the C— H bonds in propane and the alkyl 
groups, and 1.08 for the C — H bonds on the 
benzene ring. Tetrahedral angles were assumed 
in propane and in the alkyl groups. In every case 
the numerical values of the elements of the 
kinetic-energy matrices (tables 2 and 3) have the 
units of mass in atomic-weight units multiplied by 
the square of a length in angstrom units. 

Table 2. — Kinetic-energy matrices for propane before 
removal of over-all and internal rotations 



Table 3. — Kinetic-energy matrices for propane, toluene, and 
ethylbenzene 





Ax 


1 « 


<r~ a 


1.960 -2.04 
-2.04 11.63 







4>z 


0+ 


03 


7T + 




! 0* 


58. 44 


2.61 


1.60 


4.26 


A 2 


1 + 


2.61 


3.194 








I 03 


1.60 





1.597 







{ *+ 


4.26 








1.996 





Propane (C2r) 








a~ 


a 


Ax 


f -- 


1.96 


-2.04 


i « 


-2.04 


11.6 






03 


7T" 1 " 


At 


f *8 


1.55 


-0.120 


I ** 


-0. 120 


1.67 






<r + 


<ji 


Bi 


( ,♦ 


1.65 


-0.084 


I n 


-0. 084 


.724 






7T _ 


7T.3 


Bi 


[ x- 


1.63 


0.44 


[ „ 


0.44 


1.62 



Toluene (C») 


A' 


I :; 


<71 


<T2 


1.87 
-1.90 


-1.90 

14.4 


A" 


I " 

[ 7T2 


7T1 


7T2 


1.91 
-2.37 


-2.37 
20.6 



( 







0V 


«r+ 


0-3 




f *. 


66.66 


4.67 


1.57 


B, 


" + 


4.67 


1.978 


0.026 




I <T3 


1.57 


0.026 


.761 







0.r 


0- 


7T 


7T3 




f 0X 


17.80 


3.69 


-2.07 


3.16 


B 2 


0~ 


3.69 


3.194 








7T~ 


-2.07 





1.941 


-0. 045 




{ TTS 


3.16 





-0. 045 


2.358 



Ethylbenzene (C u ) 


A' 


1 <T2 

C3 

I a 


ai ff2 oz a. 


1.93 -1.85 -0.023 -2.77 
-1.85 29.8 -.57 19.9 
-0.023 -0.57 .764 -0.157 
-2.77 19.9 -.157 27.2 


A" 


03 

7T1 

7T2 

I 7T3 


03 7T1 7T2 7T3 


1.57 -0.0S9 -0.42 0.067 

-.089 1.66 .51 .25 

-.42 .51 22.4 -3.14 

.067 .25 -3.14 2.04 



I 
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As stated previously, the intermediate steps in 
i lie calculation are dependent upon the choice of 
sc,y,0-axes, but the final kinetic-energy matrices 
(after removal of over-all rotation) are not. The 
calculations for propane are described in some 
detail in order to illustrate the method, and the 
choice of axes for propane will therefore be stated. 
The x- and s-axes were taken in the plane of the 
C — C — C angle, with the origin at the vertex of 
the angle, and the :-axis bisecting the angle. The 
--axis therefore remains a twofold axis of sym- 
metry as the C — C— C angle is deformed. The 
origin is not at the center of mass of the molecule. 

The coordinate a in propane is most easily 
treated as follows. By definition the e-axis 
bisects the C — C — C angle. Denote the changes 
in these two half-angles by a x and a 2 . The 
kinetic-energy matrix is first set up for a x and a 2 , 
and then the following transformation is carried 
out : 



a=(ai + a 2 ) 
Q=(ai — a 2 ) 



(ID 



The second equation is merely the condition that 
the 2-axis bisect the C — C — C angle when the 
latter is deformed. The elements in the kinetic- 
energy matrix for (a { — a 2 ) are then ignored. The 
nonorthogonal matrix of this transformation is 



■I 

i j 



Equation 26 leads to the relations 

(a,a) = (l/4)[(a l ,a l )+2(a h a 2 ) + (a 2 ,a- 2 )] (78) 

(a,2) = (l/2)[(a 1 , 2 ) + (a 2;2 )], 

where q is any coordinate other than a h a 2 , or a. 
The following remarks apply to the calculation 
of the elements of the unreduced kinetic-energy 
matrix, eq 40, for propane. If the origin of the 
£, v, T-axes (the point x or , y or , z or ) is located at the 
carbon atom on the £-axis (or at any point on the 
£-axis for the balanced rotations <fo, </> 2 , and <fe), 
then It v and 1^, vanish in each case, and H^, 
H { y?, and H ( £ may be calculated from eq 61. 
H^, Hfr, and H^ vanish for the balanced rota- 



tions 0i, <&>, and 03. The diagonal terms K rr are 
obtained from eq OS. For balanced rotations 
K rr =I{.. Coordinates with different subscripts 
(1, 2, or 3) are unconnected, and K ry is calculated 
from eq 69. K r y vanishes if either rotation is 
balanced (^-coordinates), or if the directions of 
motion of the centers of mass are perpendicular 
(thus ^-coordinates are perpendicular to a- and 
a-coordinates) . Coordinates with the same sub- 
scripts are connected, and K ry) is obtained from 
eq 70 (the first term being given by eq 74). 
From a consideration of the directions of the 
£-axes, and the directions of motion of the centers 
of mass, for each coordinate, it is readily seen that 
the only terms that do not vanish are K ai{ai) 
tf=l,2). 

The numerical values of the elements of the 
matrices of eq 40 (one matrix for each representa- 
tion) are given for propane in table 2. 

The next step is the removal of the over-all 
rotations by the reduction process of eq 16. hi 
the subsequent calculations the frequencies of the 
0, and 0o internal coordinates (methyl rotations in 
propane and toluene, and methyl and ethyl rota- 
tions in ethylbenzene) have been assumed to be 
zero. This is justifiable because the potential 
barriers for these rotations 10 correspond to suffi- 
ciently low frequencies for small oscillations from 
the equilibrium configuration. The coordinates 
0i and 02 have, therefore, also been removed by 
the reduction process of eq 16. 

The final reduced kinetic-energy matrices, one 
for each representation, for each molecule, are 
shown in table 3. 

2. Calculation of Frequencies 

The force constants are now to be evaluated 
from the observed frequencies of propane and 
toluene. These force constants may then be sub- 
stituted in the secular equation for ethylbenzene 
to yield calculated frequencies for the latter 
molecule. 

The frequency assignment for propane has been 
taken from Pitzer [9], and the assignment for 
toluene from Pitzer and Scott [10], The fre- 
quencies (v) are summarized in table 4. The 



10 Propane, reference [9]; toluene and ethylbenzene, reference [11]. 



Frequencies of Semirigid Molecules 



13 



normal coordinates are mixtures of the coordinates 
in each representation, and when this mixing* is 
appreciable it is not possible to assign the fre- 
quencies in each representation to the coordinates 
without some ambiguity. However, the pairing 
indicated is fairly certain. Previous investigators 
have assigned the two B 2 frequencies as t~ (748) 
and 7T 3 (1179), but the assignment of these fre- 
quencies in table 4 is supported by (1) the force 
constants obtained for fa, o- 3 , and 7r 3 , which should 
have approximately the proportions 2,1, and 1, 
respectively (these figures assume that interaction 
constants are negligible, and that the main con- 
stants arise from the constants for the four 
C — C — H angles), and (2) the approximate 
equality of the force constants obtained for w~ 
and o~~. 

Table 4. — Frequencies and calculated force constants for 
propane and toluene 



PROPANE 


Coordinate 


Description 


Sym- 
metry 


V 


X 


k 


0— 


CH3 rocking . ... ... 


Ax 
Ax 
Ai 
Ai 
Bx 
Bx 
B2 
B 2 


cm-i 
1155 
375 

1,278 

940 

922 

1,338 

1,179 

748 


1.334 
0. 1406 
1.633 
0.884 
.850 
1.790 
1.390 
0.560 


2.08 
1.68 
2.50 
1.49 
1.41 
1.28 
1.97 
0.96 


C-C-C bending 

CH2 twisting . ._ 


03 

7T+" 

<r+ 

<7"3 

7T~ 

7T3 


CH3 rocking _ 


CH3 rocking 


CH2 rocking . 


CH3 rocking . ._ 


CH2 rocking 




TOLUENE 


<T\ 

<n 

7T1 

7T2 


CH3 rocking 


A' 
A' 
A" 
A" 


1,190 
216 

1,060 
340 


1.416 
0. 0467 
1.124 
0. 1156 


2.28 
0.68 
1.80 
2.43 


-CH3 rocking 8 

CH3 rocking 


-CH3 rocking 9 



The next column of table 4 gives the value of 
X (eq 22 to 24) in the units 

X-Kcm-O/lOOO] 2 . 

The last column shows the calculated values for 
the force constants, k, as calculated from the 
secular equations. They are the diagonal ele- 
ments in the potential energy, or F, matrix 
(off-diagonal elements have been set equal to 
zero as there are not sufficient data to determine 
them). The units of these force constants corres- 



pond to the units that have been used for X and 
for the kinetic-energy matrices in table 3; in the 
reverse calculation for ethylbenzene the frequen- 
cies are again obtained in cm -1 . 

Table 5 shows the force constants selected for 
the ethylbenzene molecule. Those derived from 
toluene are a 2 from a 2 , and tt 2 from 7r 2 . Those 
derived from propane are a x from (o + -{- ^~) /2 , 
Ti from (7r + + 7r~)/2, and <7 3 ,7r 3 , <£ 3 and a from a Zl 
7r 3 , 03, and a, respectively. The averaged force 
constants, (a + +a~)/2 and (t+ + ir~) /2 , are the main 
force constants for methyl rocking in propar 
The observed separation of the frequencies lur 
o- + , a~, and for 7r + , t~, in propane, is due to methyl- 
methyl interaction constants and the twofold 
symmetry. 

Table 5. — Force constants and calculated frequencies for 
ethylbenzene 



ETHYLBENZENE 


Coordinate 


Description 


Sym- 
metry 

A' 
A' 
A' 
A' 

A" 
A" 
A" 

A" 


k 


X 


V 


ax 

<72 

<T3 

a. 

<f>i 

7T1 

7T2 

7T3 


CH3 rocking 


1.75 
0.68 
1.28 
1.68 
2.50 
1.73 
2.43 
0.96 


1.067 
0.0187 
1.721 
0. 1437 
1.615 
1.144 
0. 1025 
.604 


cm— > 
1,035 

137 
1,310 

380 
1,270 
1,070 

320 

775 


C2H5 rocking _ 


CH2 rocking- 


C-C-C bending 

CH2 twisting 

CH3 rocking 

C 2H5 rocking .. 


CH2 rocking 



Substitution of these force constants in the sec- 
ular equations for ethylbenzene then yields the 
values of X, and the calculated frequencies, shown 
in the last two columns of table 5. The observed 
frequencies for ethylbenzene are discussed in the 
next section. 

IV. Vibrational Frequencies of 
Ethylbenzene 

The present calculations were undertaken as an 
aid to the assignment of the vibrational frequencies 
of ethylbenzene. * These frequencies were required 
for the calculation of the thermodynamic functions 
by statistical methods [1 1]. The frequencies have 
been assigned on a semiempirical basis which 
involved a detailed consideration of only the 



*F. G. Brickwedde, M. Moskow, and R. B. Scott [24] have also presented 
an assignment of the vibrational frequencies of ethylbenzene. 
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lowest frequencies of vibration of the benzene ring 
and also of the frequencies associated with the 
vibration of the ethyl group. The remaining 
frequencies w re taken from the assignment for 
toluene [10]. A complete frequency assignment 
for ethylbenzene at this time is impracticable. 

Several workers have reported Raman spectra 
[1 A, 15, 16, 17, 18] and infrared spectra [19, 
20, 21, 22, 23] for ethylbenzene. The lowest 
observed frequencies in the Raman spectrum of 
ethylbenzene, are, approximately, 154, 305 (weak) 
99 (weak), 452, 485, 538 (weak), 553 (weak), 
and 020 cm" 1 . The frequencies 480, 567, and 
617 cm"" 1 have been reported in the infrared spec- 
trum in this region. It will be helpful to list 
the six lowest frequencies of toluene, as assigned 
by Pitzer and Scott [10]: 11, 216 (R); 18b, 340 
(R); 16a, 405 (R) ; 16b, 467 (R) ; 6a, 521 (R) ; and 
6b, 622 (R, IR). The first number designates 
the mode of vibration (see fig. 6, of Pitzer and 
Scott), the second is the frequency in cm -1 , and 
R and IR refer to the presence of the line in the 
Raman or the infrared spectrum, respectively. 

The vibrations 11 and 18b correspond princi- 
pally to the rocking of the alkyl group perpendic- 
ular to and parallel to the plane of the benzene 
ring, respectively. The frequencies of these modes 
in ethylbenzene were calculated as 137 cm" 1 for 
11 and 320 cm -1 for 18b by means of the approxi- 
mate normal coordinate analysis of section 1 1 1 
(these 1 vibrations were denoted by a 2 and t 2 , 
respectively, in that section). The only possibility 
for the assignment of the 154 cm -1 frequency is 
therefore mode 11. The mode 18b was assigned 
the calculated frequency of 320 cm" 1 . The 
frequency observed at 305 cm -1 may easily be 
the overtone of 154 cm -1 . 

The mode 16a has a frequency of 400 cm" 1 in 
benzene and 405 cm -1 in toluene. The observed 
frequency at 399 cm -1 in ethylbenzene may 
represent 16a, or an overlap of 16a with the 
carbon skeletal bending in the ethyl group, which 
should fall at about 380 cm -1 . The toluene value, 
405 cm -1 , was retained for ethylbenzene. The 
mode 16b has the frequency 467 cm -1 in toluene; 
the nearest observed frequency in the spectrum 
of ethylbenzene, 452 cm -1 , was assigned to this 



mode. The mode 6a in toluene gives rise to a rel- 
atively intense line at 521 cm -1 . The variation of 
the frequency of the mode 6a in other monosub- 
stituted benzenes, as a function of the mass of the 
attached group [12], indicates that 6a should be 
lowered somewhat in ethylbenzene, since the 
effective mass of the ethyl group in this mode is 
considerably greater than that of a methyl group. 
Therefore, the frequency 485 cm -1 , which is the 
strongest observed line in this region, was assigned 
to 6a. The mode 6b has a frequency of 606 cm -1 in 
benzene, 622 cm" 1 in toluene, and remains essen- 
tially unchanged in frequency in other monosub- 
stituted benzene 112]. The frequency observed in 
ethylbenzene at 620 cm -1 was therefore assigned 
to 6b. 

The mode 7a, which is principally the stretch- 
ing of the C — C bond joining the alkyl group to 
the ring, has a frequency of 1,210 cm -1 in toluene. 
As a strong line appears at about 1,200 cm -1 in 
the ethylbenzene Raman spectrum, the toluene 
frequency was retained for 7a. The other C — C 
stretching mode 4 in the ethyl group was assigned 
the frequency of the strong line at 965 cm -1 in the 
Raman spectrum of ethylbenzene as this is close 
to the average of the two carbon-carbon stretch- 
ing frequencies of 868 and 1,053 cm -1 in propane 
[9]. The symmetry that splits these two com- 
ponents in propane is absent in ethylbenzene. 

The symmetrical hydrogen bending in the 
methyl group was assigned the frequency 1,385 
cm -1 , and the doubly degenerate unsymmetrical 
bending the frequency 1,460 cm -1 . The symmet- 
rical hydrogen bending in the methylene group 
was assigned the value 1,460 cm -1 . The value 
2,960 cm -1 was assigned the C — H stretchings in 
the methyl and methylene groups. These assign- 
ments were made on the basis of the frequencies of 
these modes in propane and toluene, and to some 
extent on the basis of the spectrum of ethylben- 
zene. 

The remaining frequencies of vibration of the 
ethyl group were assigned values calculated by 
means of the approximate normal coordinate 
analysis of section III (see tables 1 and 5) . 

The complete frequency assignment for ethyl- 
benzene is shown in table 6. 
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Table 6. — Vibration frequencies of the ethylbenzene molecule 



Type of Vibration 



Sym- 
metry 



Fre- 
quency 



Vibrations principally characteristic of the benzene ring or phenyl radical 



11 (C 2 H 5 rocking).. 
18b (C2H5 rocking) _ 

16a 

16b 



6b— 
4._.. 
10b- 
12_._ 
10a. . 
17b- 
5.... 
17a_. 
1.... 



15- 
9b. 



7a (C6H5-C2H5 stretching). 



19b.. 
19a_. 
8b— 
8a— 
14... 
13 

7b 

2 
20a 
20b 



>C-H stretching.. 





cm 1 


A' 


154 


A" 


320 


A" 


405 


A' 


452 


A' 


485 


A" 


620 


A' 


695 


A' 


730 


A' 


785 


A" 


842 


A' 


890 


A' 


943 


A" 


985 


A' 


1,002 


A' 


1,030 


A" 


1,070 


A" 


1,155 


A' 


1,175 


A' 


1,210 


A" 


1,282 


A" 


1,310 


A' 


1,483 


A" 


1,586 


A' 


1,603 


A" 


1, 630 


A' 


3,046 


A" 


3,047 


A' 


3,062 


A' 


3,080 


A" 


3,080 



Vibrations principally characteristic of the ethyl group 



C6H5-CH2-CH3 angle bending. 

CH2 rocking 

CH2-CH3 stretching 

CH3 rocking 

CH3 rocking 

CH2 twisting 

CH2 rocking 

CH3 symmetrical bending 

CH2 symmetrical bending 

CH3 unsymmetrical bending.. 

CH2 symmetrical stretching.. _ 
CH 2 unsymmetrical stretching 
CH3 symmetrical stretching _ _ 

CH3 unsymmetrical stretching 



A' 


380 


A" 


775 


A' 


965 


A' 


1,035 


A" 


1,070 


A" 


1,270 


A' 


1,310 


A' 


1,385 


A' 


1,460 


A' 


1,460 


A" 


1,460 


A' 


2,960 


A" 


2,960 


A' 


2,960 


A' 


2,960 


A" 


2,960 
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Appendix 1 n 

Let the coordinates q\ be divided into two groups, q a 
and qb, and partition Sn into the corresponding subma- 
trices. Then the kinetic energy matrix, eq 4, becomes 

Saa 8ab S a -2 



S a b 



Sbb 
Sh2 



Sv, 



(79) 



_S a 2 Sb2 $22_ 

The removal of the cross terms in the kinetic energy that 
contain q a leads to the reduced matrix 

(S bb — S^Sa^Sai) ($ b2 — ^ab^al 1 ^) 

_ (£&2 — S a2 S a ~ l S ab ) ($22 — S^Sa^S^) _ 

and the removal of the remaining cross terms that contain 
qb leads to the final reduced matrix for the coordinates #2, 

11 See section II, 1 (following eq. 16). 



D 



;] 



(80) 
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(S b2 -^8- l S a2 ). 
Equation 81 may be rewritten in the form 

L — G bb S ab S a ~ l 

(*bb J LA2 J 



(81) 



~£>aa 1 $abGbb 



where 



Gbb — (Sbb ~ ^ab^aa 1 S ab ) 



(82) 



(83) 



On the other hand, the simultaneous removal of the cross 
terms that involve q a and q h leads to the reduced matrix 



S 22 =S 22 -[S t 



U H* S 4'M (84) 

LS ab s bb J Lv>_ 



The matrices $' 22 calculated from eq82and 84 are identical, 
provided that 



(S -1 + £ -»£ h Gu h S h S - 1 ) — ^ 
— P'S! .Sf _1 



Cm J 



["&„ S«f| r^aa "I 

LS«* S«J LO R J 



(85) 



where E aa and Em arc unit matrices. Bui the last rela- 
tion is easily verified by multiplying out the matrices on 
the left. This establishes the theorem when r/, is divided 
into i wo groups; the truth of the t heorem for any number 
of groups follows by induction. 



Appendix 2 l2 



From the relation that SG = E, where K is the unit 
matrix, and eq 3, there are obtained the relations 



S 12 G l2 + £»Gn + £ 2 3&23=E 
5fll(?12 + ^ 12 (? 22 + ^3(?23 = 

S\ 2 G\z -h S 22 G 2 3 -\- SmGm = 
SnGiz + Si 2 G 2 z + «Si3(j33 = 



(86) 

(87) 

(88) 
(89) 



Now multiply eq 89 from ihe left by SumSiT 1 , and from the 
right by G& G 23f and add t lu^ resulting equation to eq 
86; next multiply eq 87 from the left by SwSu 1 and eq 88 
from the right by Crsi" 1 ^, and subtract both of the result 
ing equations from the previous sum. The resulting 
equation may be factored as follows: 

(S 22 -S i2 Sn l S 12 ) (G 22 - G u G» l Qj = E 22 . (90) 

Equation 21 follows immediately from eq 16, 20, and 90. 

Washington, September 24, 1946. 



u See section II, 1 (eq 20 and 21) 
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